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ORIGINAL ARTICLE 



A gorilla reservoir for human T-lymphotropic virus type 4 

Matthew LeBreton^'^'^'*, William M Switzer^'*, Cyrille F Djoko^, Amethyst Gillis^'^, Hongwei Jia*, 
Michele M Sturgeon^, Anupama Shankar^, Haoqiang Zheng*, Gerard Nkeunen^, Ubald Tamoufe^'^, 
Ahmadou Nana^, Joseph Le Doux Diffo^, Babila Tafon^, John Kiyang^, Bradley S Schneider^, Donald S Burke'^ 
and Nathan D Wolfe^'*'^ 

Of the seven known species of human retroviruses only one, human T-cell lymphotropic virus type 4 (HTLV-4), lacks a known animal 
reservoir. We report the largest screening for simian T-cell lymphotropic virus (STLV-4) to date in a wide range of captive and wild 
non-human primate (NHP) species from Cameroon. Among the 681 wild and 426 captive NHPs examined, we detected STLV-4 infection 
only among gorillas by using HTLV-4-specific quantitative polymerase chain reaction. The large number of samples analyzed, the diversity 
of NHP species examined, the geographic distribution of infected animals relative to the known HTLV-4 case, as well as detailed 
phylogenetic analyses on partial and full genomes, indicate that STLV-4 is endemic to gorillas, and that rather than being an ancient virus 
among humans, HTLV-4 emerged from a gorilla reservoir, likely through the hunting and butchering of wild gorillas. Our findings shed 
further light on the importance of gorillas as keystone reservoirs for the evolution and emergence of human infectious diseases and provide 
a clear course for preventing HTLV-4 emergence through management of human contact with wild gorillas, the development of improved 
assays for HTLV-4/STLV-4 detection and the ongoing monitoring of STLV-4 among gorillas and for HTLV-4 zoonosis among individuals 
exposed to gorilla populations. 
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INTRODUCTION 

Human infection with a plethora of simian retroviruses is well docu- 
mented and has led to global pandemics, exemplified by the human 
immunodeficiency virus (HIV) and human T-lymphotropic virus 
(HTLV), and local outbreaks, as occurs with simian foamy virus 
(SFV), simian T-lymphotropic virus (STLV) and occasionally simian 
immunodeficiency virus (SIV).'"' Given the demonstrated pandemic 
potential of retroviruses, a full understanding of the epidemiology and 
animal reservoirs of zoonotic primate retroviruses is of importance for 
monitoring and preventing future retrovirus pandemics. 

The primate origin of these zoonotic infections has been identified 
by detailed epidemiological and phylogenetic analyses of both human 
and non-human primate (NHP) retroviruses. Through detailed exam- 
ination of SIV infections in over 45 NHP species, it has been demon- 
strated that HIV-1 arose from multiple cross-species transmissions of 
SIVcpz and SIVgor from chimpanzees and gorillas, to people in west 
Central Africa. ^'^ Similarly, HIV-2 in West Africa originated from 
multiple introductions of related SIVs from sooty mangabeys.'"^ 
Many studies have described ongoing zoonotic infection of primate 
workers, hunters and butchers with SFV from a variety of monkeys 
and apes and phylogenetic analysis showing that the co-evolution 
of SFV with NHPs has facilitated accurate identification of the 
simian origin of infection.^"'" Likewise, phylogenetic analysis has 



demonstrated that at least three of the four major HTLV lineages, 
HTLV-1, HTLV-2 and HTLV-3, originated from multiple introduc- 
tions of their simian virus counterparts in monkeys and apes, STLV- 1 , 
STLV-2, and STLV-3, respectively.''*"** Thus, STLV and HTLV 
lineages are called primate T-lymphotropic viruses (PTLVs). 

Of the four major HTLV lineages, the majority of human infec- 
tions are caused by HTLV-1, which is conservatively estimated to 
be responsible for 5-10 million global human infections." HTLV- 
2, while less frequent, has also spread globally and like HTLV-1, is 
transmitted from mother to child through breast feeding, sexual 
contact and contaminated blood products during transfusion or 
injection drug use." HTLV-1 causes adult T-cell leukemia, HTLV- 
associated myelopathy/tropical spastic paraperesis and other 
inflammatory diseases in 2%-5% of infected persons. '^'^' HTLV- 
2 is less pathogenic, but is increasingly associated with HTLV- 
associated myelopathy/tropical spastic paraperesis and some 
inflammatory diseases.^* HTLV-3 is much less frequent and has 
only been found in people in Cameroon who have direct contact 
with NHPs from which they were likely infected.*"**''^'"^ No direct 
evidence of disease or secondary transmission of HTLV-3 has been 
reported.^'*"* 

The fourth HTLV lineage, HTLV-4, was discovered in 2005 in a 48- 
year-old male hunter (1863LE) from rural Cameroon who reported 
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hunting monkeys, chimpanzees, gorillas and other animals. '^^ The 
infection of individual 1863LE currently represents the only known 
human or simian virus in this lineage. As for HTLV-3, person-to- 
person transmission and disease have yet to be reported for HTLV-4, 
though both viruses possess genes and motifs believed to be associated 
with disease potential in HTLV-1 and HTLV-2.'''''^'** Molecular dating 
suggests that HTLV-4 is the oldest PTLV lineage, having originated 
some 200 000 years ago (ya).'^ The ancient age of the HTLV-4 lineage 
led to the hypothesis that HTLV-4 is a descendent of an ancestral PTLV 
that infected humans during their evolutionary history and frequent 
NHP contact in Africa and represents a rare strain circulating within 
people in the central African region.'^ Alternatively, it has been 
hypothesized that HTLV-4 originated from a more recent zoonotic 
infection with a divergent STLV also present in NHPs cohabitating 
the forests of Cameroon.^* Identification of the simian origin of 
HTLV-4 would help determine the validity of each hypothesis. 

WhUe STLV has been found in at least 30 African and Asian NHPs, 
discovery of the simian counterpart of the recently identified divergent 
HTLV-4 lineage found in a hunter from Cameroon has been elusive, 
despite surveys of potential monkey reservoirs.'^"^' In an attempt to 
discover how HTLV-4 emerged, and to shed light on its evolution and 
pandemic potential, we screened a large number of monkey and ape 
blood specimens from Cameroon using a recently developed HTLV-4- 
specific quantitative polymerase chain reaction (qPCR) assay.^^ 

MATERIALS AND METHODS 

Ethics statement 

Specimen collection from captive and wild primates was approved by 
the University of Johns Hopkins Animal Care and Use Committee 
(FS03M221 and FS06H205), University of Cahfornia Institutional 
Animal Care and Use Committees (lACUC) (ARC#2007-110-01A, 
ARC#2007- 110-2, ARC#2007-110-3), the University of California 
Davis lACUC (Protocol #16067) and the Cameroon government. 
All sample collections from captive primates were undertaken by ve- 
terinary animal care staff (authors Babila Tafon and John Kiyang) in 
the Limbe WQdlife Centre and Mfou National Park (wUdlife rescue 
centers) during immobilization as part of quarantine health checks, 
regular health checks or during transport between enclosures. 

Specimens from hunted primates were obtained by opportunistic 
sampling of animals in the wildlife trade or animals killed by local 
hunters as part of an ongoing effort to gather samples from bushmeat 
for detection of viruses to help elucidate the risks of this activity. 
Generally, these samples were self-collected by the hunter and gath- 
ered by members of our team at a later time point. Hunters never 
received money for specimens and particular care was taken to fully 
avoid promoting these activities. Our team had absolutely no control 
over the species selected by hunters, and also worked closely in the field 
with staff from the Cameroonian ministry responsible for wildlife to 
ensure that it was clear to hunters which species were protected and 
should not be hunted. Nonetheless, a small number of protected spe- 
cies were hunted in these sites; these samples were obtained and tested 
along with other samples. We work with the communities in 
Cameroon and regionally in Central Africa to address such issues 
and actively discourage hunting of protected species in the sites where 
we work. All protocols were developed in partnership with and with 
guidance from leading conservation organizations, including the 
Wildlife Conservation Society, the University of California Davis 
and EcoHealth Alliance (formerly, the Wildlife Trust) to ensure the 
highest ethical standards as we conduct trainings and build wildlife 
surveillance capacity. 



Sample collection and preparation 

Specimen collection from captive and wild animals was approved by 
lACUC and the Cameroon government. People who reported hunting 
or butchering wild animals in rural areas in Cameroon volunteered to 
collect blood samples from wild animals that they had hunted.^" 
Hunters were educated about the risks associated with direct contact 
with wUd animals and about appropriate prevention measures. Dried 
blood spots were collected on Whatman paper during butchering and 
allowed to air dry before storage in labeled envelopes. Preliminary 
identification of hunted species was undertaken by using pictures of 
animals common in the region and lists of local names of various 
hunted species were compiled.^' Host species identity was confirmed 
for positive specimens using analysis of mitochondrial sequences 
(detailed below). Additional blood samples were obtained by techni- 
cians directly from hunted animals when encountered in the rural 
villages, as well as from captive wild animals orphaned by the bush- 
meat trade and housed at two sanctuaries in Cameroon. Samples from 
captive animals were obtained during routine veterinary procedures 
such as quarantine health checks, recaptures and during illness. Whole 
blood samples obtained by technicians from captive animals were 
collected in EDTA tubes and processed into plasma and buffy coat 
within 48 h and frozen at ^ — 20 °C. AU dried blood spot samples were 
stored in clip lock plastic bags with silica gel at room temperature since 
collection. Nucleic acids were extracted manually using the Qiagen 
QIAamp DNA minikit or with the Invitrogen Dynabead SILANE viral 
nucleic acid kit according to manufacturer's protocols. DNA integrity 
was confirmed by P-actin PCR as described elsewhere.^' All animal 
samples were stored in laboratory spaces where only animal specimens 
were kept. All DNA preparation was conducted in spaces dedicated to 
animal samples. DNA concentrations were measured using a 
Nanodrop spectrophotometer. 

PCR detection and generation of complete viral genomes 

In total, 5 |iL of each extracted DNA (~100 ng) specimen was first 
screened using a specific qPCR assay created to detect and quantify 
HTLV-4 polymerase (pol) sequences in a specimen.^^ Positive samples 
were confirmed using nested PCR assays employing primers based 
on the HTLV-4 genome (GenBank accession number EF488483);'^'^'' 
primers 1863PF5 (5'-ACC ACT ACC CTT GCC TCG CCC TCC-3') 
and 1863PR1 (5'-AAT GGG GAT GGT GGG GGT GGT TTC-3') and 
1863PF1 (5'-TCC AGA TTC CCC TCG CAA AGG GAT-3') and 
1863PR2 (5'-ACG ATG AGG GAT GTG GGG AAG GCT-3') amplify 
488-bp and 195-bp po/ sequences, respectively. In addition, fax sequences 
were detected using the nested PCR primers 1863TF1 (5'-CTC CTT CTT 
TCA GTG GGT GCG GAG-3') and 1863TR1 (5'-GGG GTA GTG AGG 
TTT GGG TGG TAT-3') and the internal primers 1863TF2 (5'-CCT 
ACC GCA AGG GAT GTG TTG AAA-3') and 1863TR2 (5'-TAT GGG 
GCC GGT GTG ATG ATA AAG-3') to generate 399-bp and 275-bp 
fi'agments, respectively. Partial long terminal repeat (LTR) sequences 
524-bp in size were amplified using the nested PCR primers 1863LF1 
(5'-GCG AGA GCA AGG CCC CAA GGA CAA-3') and 1863LR1 (5'- 
GGA GGG GTT TGA GTA CAG GGG GCT-3') in the first round and 
1863LF4 (5'-ACC TCG CCC TTA CCC ACT TCG CCT-3') and 
1863R1 in the second round of PCR. Nested PCR was also used to 
generate 319-bp envelope (env) sequences using generic PTLV primers 
FGENVF2 and PGENVR2 described previously.'^ Standard PCR con- 
ditions were used for 35 cycles in each round for all PCR assays. 
Complete STLV-4 genomes were obtained using high fidelity Taq poly- 
merases and long PCR with primers and conditions previously used to 
obtain the fiiU-length genome of HTLV-4.'^ 
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NHP species confirmation, gorilla mitochondrial haplotype 
determination and sex determination 

NHP identity was confirmed for positive individuals by amplification 
and sequencing of the mitochondrial DNA (mtDNA) cytochrome 
oxidase subunit II gene that was previously shown to be an excellent 
marker for primate phylogeny."* GoriUa population affinity was deter- 
mined by phylogenetic analysis of mtDNA hypervariable region 1 
(HVl) sequences and comparison to documented population vari- 
ability.^' New HVl sequences and those available from gorillas at 
GenBank were aligned with ClustalW in MEGA5, edited manually, 
and the best model of nucleotide substitution determined to be 
Tamura-Nei and gamma distributed rates. Neighbor-joining trees 
were inferred using MEGA5 with tree topology support determined 
using 1000 bootstrap replicates. Nuclear mtDNA HVl sequences were 
determined using this initial analysis and were then excluded from the 
final analysis. Sex of the wild gorillas was determined using a com- 
bination of three PGR primers in a single reaction (UTXY (forward), 
5'-TGC TAG CTG AGG TGG ACA AGA AGG-3', UTY (reverse), 5'- 
TGG TTG TTT GAG GGA GCA AGG RTG TAT K-3' and UTX 
(reverse), 5'-GTG GAG AGT GGG AGT GCT GTT AGG-3') to amp- 
lify the ubiquitously transcribed tetratricopeptide repeat protein gene 
as described in detail elsewhere.^* The assay detects tetratricopeptide 
repeat protein sequences in humans and NHPs and produces two 
fragments of 127-bp and 86-bp in size for males and a single 127-bp 
sequence for females. 

PTLV phylogenies 

Nucleotide aUgnments were made using GlustalW in the alignment 
editor program in MEGA5 and edited manually. Percent nucleotide 
identities were calculated using alignments of selected sequences using 
the software Geneious Pro v6.0.3. The best fitting distance model of 
nucleotide substitution for each alignment was inferred using the max- 
imum likelihood method with goodness of fit measured by the Bayesian 
information criterion in MEGA5. The best fitting nucleotide substi- 
tution model for the phylogenetic alignments was inferred to be the 
Hasegawa-Kishino-Yano model with discrete gamma and invariant 
among-site rate variation. Phylogenetic relationships were inferred using 
Bayesian analysis with the BEAST vl.6.2^^ program. Statistical support 
for the inferred Bayesian trees was assessed by posterior probabilities. For 
the Bayesian analyses, a relaxed or strict molecular clock model was used 
and each run consisted of two independent 100X10^ Markov chain 
Monte Garlo generations with sampling every 10000 generations and 
a constant or Yule coalescent tree prior. Gonvergence of the Markov 
chain Monte Garlo was assessed by calculating the effective sampling size 
of the runs using the program Tracer vl.5 (http://beast.bio.ed.ac.uk/ 
Tracer). All parameter estimates showed significant effective sampling 
sizes >200. The tree with the maximum product of the posterior clade 
probabilities (maximum clade credibility tree) was chosen from the 
posterior distribution of 10001 sampled trees after burning in the first 
1000 sampled trees with the program TreeAnnotator version 1.6.2. The 
most recent common ancestor dates (TMRCAs) and nucleotide substi- 
tution rates were co-inferred using BEAST and a strong prior of the 
divergence time of 40 000-60 000 ya for the Melanesian HTLV-1 hneage 
(HTLV-lmel) and 15 000-30000 ya for the most recent common 
ancestor of HTLV-2a/HTLV-2b native American strains as previously 
described.'"'^^'^'"^'' Such diver gence dates were based on well-established 
genetic and archaeological evidence, suggesting that ancestors of indi- 
genous Melanesians and Australians migrated from Southeast Asia or 
the introduction of ancestral indigenous Indians into North America via 
the Bering Strait during those times.'^"'^'^^^^ The use of two calibration 



points has previously been shown to provide more reliable estimates of 
PTLV substitution rates and divergence dates than a single calibration 
date.'^'"*'^^^** Gomparisons of the tree likelihoods for each tree prior and 
clock model were performed using Bayes factors calculated in BEAST. 

Serology 

Gross-reactive HTLV antibodies were detected in plasma samples 
available from PGR-positive captive gorillas using the HTLV Blot 
2.4 (MP Diagnostics, Santa Ana, GA, USA) Western blot kit following 
the manufacturer's instructions. The HTLV Blot 2.4 kit contains as 
antigens an HTLV-1 crude lysate, a recombinant envelope protein 
(Env, GD21) that most PTLVs cross react with and type-specific pep- 
tides that differentiate between HTLV-1 and HTLV-2. This kit has 
been used to determine the Western blot reactivity of plasma from 
the HTLV-4-infected person 1863LE'^''" and plasma from HTLV- 
4(1863LE) was also run for comparison. 

RESULTS 

Prevalence and geographic range of STLV-4 infection in wild and 
captive gorillas from Cameroon 

To broadly screen simians for potential STLV-4 infection, samples 
from a total of 1107 individual NHPs (681 wild and 426 captive) 
consisting of 21 species were screened for evidence of PTLV-4 infec- 
tion using a qPCR assay specific for the detection of HTLV-4-like pol 
sequences (Table 1).^^ The pol qPGR assay has been previously shown 
to have a dynamic linear range of nine logs and can reliably detect at 
least 10 copies of HTLV-4 in 1 ug DNA.^^ Of the 21 species tested, 
DNA samples from only three wild (09584, 09948 and 19200) and 
three captive (Benito, Nyum and Shai) gorillas were positive using 
this assay (Table 1). Thus, the prevalences of STLV-4 in captive and 
wild gorillas in our study population were 8.3% and 13.6%, respec- 
tively. BLAST analysis of mtDNA cytochrome oxidase subunit II 
sequences confirmed the gorilla {Gorilla gorilla, ggo) origin of all six 
PTLV-4-positive animals. All positive wild gorillas were found to be 
males using sex-specific PGR, while the sex of the three positive captive 
gorillas was confirmed using the same assay with only Nyum being 
female (Table 2). 

Figure 1 shows the sampling sites for the three wild STLV-4-positive 
gorillas, the origin of the three captive gorillas and the sampling lo- 
calities of the STLV-4-negative NHPs. All sbc STLV-4-positive ani- 
mals originated in East Gameroon in proximity to the village of the 
HTLV-4-infected individual. The positive samples from hunted 
gorillas originated in villages 12 km, 26 km and 55 km from the 
village of the HTLV-4-infected individual outside Lomie. The exact 
origins of the captive animals are unknown as they may have been 
moved from the area where they were originally captured. Nyum 
was confiscated from an illegal ape dealer in Belabo (235 km north 
of person 1863LE), Shai was confiscated from a hunter in Abong 
Mbang (150 km north of person 1863LE) and Benito came from a 
forestry camp near Mindourou (90 km northeast of person 1863LE) 
(Figure 1). Analysis of mtDNA HVl sequences showed that the 
origin of all six gorillas were within the range of haplotypes pre- 
viously observed in Gameroon and whose diversity is restricted by major 
rivers within Gameroon (Figure 2).^^'^''' For example, Nyum is haplotype 
Gl which occurs north of the Dja and Sanaga Rivers; Shai, Benito and 
wild gorilla 19200 are all haplotype G2 which has been shown to exist 
between the upper Dja River and the Sanaga River, while wild gorillas 
9584 and 9948 are haplotype D2 which occurs south of the Sanaga River 
in Gameroon and near the Gameroon/Central African Republic/ 
Republic of Gongo border (Figure 1)}^'^^ These results are consistent 
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Table 1 Prevalence of STLV-4 in NHPs from Cameroon 



WildNHPs Captive NHPs 







STLV-4 positive/tested 
(prevalence) 


STLV-4 positive/tested 
(prevalence) 


Hominidae 








Gorilla gorilla 


Western gorilla 


6/ZZ {Ld.b/o) 


3/3o [G.d /o) 


Pan troglodytes 


Chimpanzee 


0/24 


0/147 


Colobidae 








Colobus guereza 


Guereza colobus 


0/16 


0/1 


Colobus satanas 


Black colobus 


0/42 




Piliocolobus preussi 


Preuss's red colobus 


0/5 




Cercopithecidae 








Cercocebus agilis 


Agile mangabey 


0/22 


0/18 


Cercocebus torquatus 


Red-capped monkey 


0/46 


0/4 


Cercopithecus cephus 


Moustached monkey 


0/76 


0/5 


Cercopithecus erythrotis 


Red-eared monkey 


0/123 


0// 


Cercopithecus Fhoesti preussi 


Preuss's monkey 


0/4 


0/7 


Cercopithecus mona 


Mona monkey 




0/17 


Cercopithecus neglectus 


DeBrazza's monkey 


0/4 


0/2 


Cercopithecus nictitans 


Spot-nosed monkey 


0/141 


0/17 


Cercopithecus pogonias 


Crowned monkey 


0/35 


0/1 


Cercopithecus tantalus 


Tantalus monkey 




0/24 


Erythrocebus patas 


Patas monkey 




0/5 


Lophocebus albigena 


Grey-cheeked mangabey 


0/40 


0/3 


Mioplthecus ogouensis 


Northern talopoin 


0/1 


0/1 


Mandrillus leucophaeus 


Drill 


0/22 


0/49 


Mandrillus sphinx 


Mandrill 


0/58 


0/27 


Papio anubis 


Baboon 




0/55 




Totals 


3/681 (0.4%) 


3/426 (0.7%) 



with the presumed origin of the six positive gorillas and also verify the 
specimens are from six different individuals. 

Epidemiological and clinical considerations of STLV-4 in captive 
gorillas 

Nyum and Shai reside in the same enclosure in the same sanctuary near 
Yaounde, while Benito is in a sanctuary in Limbe and has had no known 



contact with the other two animals, or any individuals from their group. 
Shai arrived at the sanctuary in 2002; a blood sample from 2009 was 
negative by qPCR testing and subsequent samples at three separate time 
points in 201 1 were qPCR-positive (Table 2). Nyum arrived in April 2003 
and has shared the same enclosure with Shai since September 2003. 
Nyum tested qPCR-positive in 2009. No earlier samples are available 
for either of these animals. Nyum's first observed sexual encounter was 



Table 2 Identification of PTLV-4 in gorillas from Cameroon 



Sample 


Origin^ 


Sex 


Birth date 


Date collected 


Western blot profile 


PTLV-4 copies/100 
cells" 


nPCR pol, env, TL, LTR, 
FLG" 


9584 


WM 


Male 


UnknovKn 


30 Dec 2004 


NA"^ 


BLD"^ 




+ , +, +, 


4- 


9948 


Wild 


Male 


Unknown 


07 Apr 2005 


NA 


6439 




+ , +, +, 


4- 


19200 


Wild 


Male 


Unknown 


03 Nov 2007 


NA 


58 




+ , +, +, 


4- 


50539 (Benito) 


Captive: LWC 


Male 


-Apr 1996 


24Jun 2009 


HTLV-positive, 


56 




+ , +, +, 


4- 












untypeable 










51323 (Nyum) 


Captive: MNP 


Female 


-Oct 2001 


21Feb2009 


Indeterminate 


3664 




+ , +, +, 


4- 


51461 (Nyum) 


Captive: MNP 


Female 




24 Jun 2010 


Indeterminate 


5243 




+ , +, +, 


4- 


51321 (Shai) 


Captive: MNP 


Male 


-Oct 2001 


21 Feb 2009 


Indeterminate 


NT 




NT' 




51573 (Shai) 


Captive: MNP 


Male 




12 Apr 2011 


HTLV-2-like 


4536 




+ , +, +, 


4- 


51647 (Shai) 


Captive: MNP 


Male 




03 Nov 2011 


HTLV-2-like 


2340 




+ , +, +, 




51657 (Shai) 


Captive: MNP 


Male 




09 Dec 2011 


HTLV-2-like 


7348 




+ , +, +, 


+ 


HTLV-4 (1863LE)^ 


Lomie 


Male 


-1954 


14 Nov 2006 


HTLV-2-like 


0.295 




NT 





= MNP, Mfou National Park; LWC, Limbe Wildlife Center. 

Average of two experiments except for 51647 and 51657 which were tested once; assumes 100 cells/ng DNA. 
■^nPCR, nested PCR; po/ (662 bp), polymerase; eni/(-319 bp), envelope; TL(227 bp), tax-UR (long terminal repeat) fragment; FLG (8793-8794 bp), full-length genome. 

NA, specimens not available for testing. 
" BLD, below the limit of detection of the assay. 
'NT, not tested. 

* HTLV-4 1863LE complete genome, GenBank accession number EF488483. 
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★ HTLV-4 1863LE '^'^^""^ 
▲ Location of confiscation of captive STLV-4 positive gorilla; 

# STLV-4 positive wild gorilla 
# STLV-4 negative non-hum 

primates 

^— IVIajor rivers 
^— IVIain roads 
— Secondary road 




Figure 1 Location of STLV-4-positive gorillas and STLV-4-negative non-human primates in Cameroon. Mitochondrial hypervariable region 1 haplotypesare included 
following animal names or specimen numbers (i.e. CI, C2 and D2). Triangles show location where captive animals were confiscated which may not be the same 
location the animals were caught in the wild. 



in June 2012; however, no activity with Shai has been observed. 
Nonetheless, these results suggest that Shai may have infected Nyum in 
captivity. Benito arrived in 1996 and a single blood sample from 2009 was 
qPCR-positive. Blood specimens from 32 additional wild-born and one 
captive-born captive gorillas, including all wUd-born animals that share 
enclosures with the positive animals at the two sanctuaries, and including 
the animals that arrived at the same time as positive individuals, were 
tested and were all found to be negative for STLV-4 sequences. 

To further investigate infection in the STLV-4-positive gorillas, we 
tested plasma samples available from the three captive gorillas using an 
HTLV-1/2 Western blot kit previously shown to detect antibodies in 
the HTLV-4-infected hunter. We found a broad range of cross-reac- 
tivity in the STLV-4-positive samples to the HTLV-1/2 antigens (Figure 3 
and Table 2). Samples from Shai and Benito were seropositive showing 
reactivity to the p24 Gag and GD21 Env proteins (Figure 3). The 2009 
sample from Shai was weakly seroreactive to the pl9 Gag protein only, 
whereas his subsequent samples in 2011 were all HTLV-2-like with 
reactivity to the HTLV-2-type specific peptide (rgp46-II) reflective of 
possible seroconversion near the first time point in 2009 with infection 
likely obtained from Nyum in captivity. Plasma from Benito showed a 
seropositive but untypeable pattern without reactivity to the rgp46 
type-specific peptides. Both plasma samples from Nyum were seroin- 
determinate and were weakly seroreactive to only the GD21 recombi- 
nant Env protein. Plasma from HTLV-4 (1863LE) was also HTLV-2-like 



as previously reported, ' but in contrast to the STLV-4-infected gori- 
llas, also showed reactivity to the p26-p32 Gag proteins. 

We also determined proviral loads (pVLs) in the STLV-4-infected 
gorillas and HTLV-4-infected hunter using our qPCR test to evaluate 
the potential role of increased pVLs in transmission and disease pro- 
gression as occurs with other retroviruses, like HTLV and HLV. pVLs 
detectable in the gorillas ranged from 56 copies/ 100 cells to 7348 
copies/100 cells (mean=3647 copies/100 cells; median=3849 copies/ 
100 cells) (Table 2). pVLs in wild-caught goriUa 9584 were below the 
limit of detection. In contrast, pVLs in the HTLV-4-infected hunter 
(1863LE) were very low (0.295 copies/100 cells) compared to the 
gorilla pVLs. The two gorillas, Nyum and Shai, which shared the same 
enclosure had comparable pVLs averaged from their multiple time- 
points (4453 and 3229 copies/100 cells, respectively), but a trend 
towards increased pVLs over time was observed in both animals. 

The potential for PTLV-4 to cause disease like HTLV- 1 and HTLV- 
2 is not clear. The only HTLV-4-infected person was apparently 
healthy in 2006; however, his family reported in 2009 that he had 
passed away after a short illness, characterized by fatigue and respira- 
tory symptoms. Animal records showed that aU three captive gorillas 
with STLV-4 infection experienced illnesses common in wild and 
captive animals, including wounds from aggressive behavior, skin 
infections, colds, fever and conjunctivitis. However, Nyum has 
developed some lameness and possible brain damage after an acute 
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Figure 2 Gorilla mitochondrial hypervariable region 1 haplotype determination. Haplotypes inferred by neigtibor-joining analysis usingan alignment of eigtit new sequences 
from ttie current study and 158 reference sequences available at GenBank. Final alignment length was 210-bp. Bootstrap support >70 is provided at branch nodes. 
Nuclear mtDNA sequences were not included in the analysis. Gorilla sequences 51461 and 51323 are from different time points from the same animal (Nyum), whereas 
the two sequences from gorilla 9948 are different clones prepared at different times to verify the robustness of the cloning procedure and phylogenetic analysis. 



illness with fatigue, fever, seizures, lethargy, deafness, loss of vision and 
temporary paralysis in 2005. It is unclear if the high pVLs in this 
animal are associated with these possible neurological symptoms. 
Blood samples were not available from 2005 for qPCR testing when 
these symptoms were reported. 

Phylogenetic analysis of STLV-4 and dating of the evolutionary 
history of PTLV-4 

Partial PTLV-4 pol, em and tax sequences were detected by nested 
PGR for the seven gorilla samples (Table 2) and all showed at least 
99% nucleotide identity to HTLV-4 (1863LE) and to each other with 
the exception of gorilla 9584 which had a PTLV-l-like env sequence 
(see below). PTLV-4 LTR sequences from all six animals, which are 
typically one of the most divergent regions in PTLV since it is a non- 
coding region and is not under the same evolutionary constraints as 
coding regions, shared 98.4%-99.5% nucleotide identity with the 
HTLV-4 LTR. The greatest LTR sequence divergence was seen within 



the wild gorillas and ranged within 97.9%-99.5% compared to 
99.8%-100% in the captive gorillas. Combined, these results confirm 
the identification of STLV-4 and demonstrate limited diversity of 
STLV-4 in different gorilla populations in Cameroon. These results 
were confirmed by phylogenetic analysis of tax and pol sequences 
which showed a clustering of all STLV-4 sequences with those from 
HTLV-4 (Figure 4). 

To further demonstrate the genetic relatedness of STLV-4 and 
HTLV-4 and to investigate changes that may occur following cross- 
species transmission, we obtained complete genomes from all three 
captive animals (Benito, Nyum and Shai), including complete ge- 
nomes from Nyum and Shai at multiple time points allowing an inves- 
tigation of longitudinal virus evolution (Table 2). All complete STLV-4 
genomes were 8793-8794 bp in length, were nearly identical to each 
other and shared 99.4%-99.7% nucleotide identity to the 8791-bp 
genome of HTLV-4 (1863LE) due to 28-50 point mutations through- 
out the genome, but not affecting the open reading frames. 1919-bp of 
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Figure 3 HTLV Western blot serologic pattern of STLV-4-infected captive gorillas 
from Cameroon. Asterisks indicate specimens witti indeterminate Western blot 
results. Positive controls are those supplied w/ith manufacturer's kit. Reactivity to 
HTLV-specific proteins is indicated on the left. 

the env-tax region was obtained from one wild gorilla (9948), which 
showed at least 99.3%-99.5% nucleotide identity to STLV-4 from the 
captive gorillas and 99.5% identity to HTLV-4 (1863LE). Five of the 
nucleotide substitutions in the complete genomes were G>A mutations 
likely introduced via ABOPEC deamination, but which did not result in 
the introduction of premature stop codons (Figure 5). Twice as many 
G>A mutations were also present in the sequences of two hunters 
found to be infected with STLV-3 most likely from a red-capped man- 
gabey (Figure 5). The HTLV-3 genome of one these individuals (Pyl43) 
also contained a 366-bp deletion in the pX proximal region that may 
also be further evidence of host restriction, but sequencing of STLV-3 
genomes from additional mangabeys is needed to verify this hypothesis. 
Nonetheless, these results suggest the absence of gross genetic muta- 
tions that would affect infectivity and replication of PTLV-4. 

To estimate the possible date of STLV-4 entry into humans, detailed 
Bayesian phylogenetic analyses were performed. Sequence analysis 
showed unequal base composition for some lineages and substitution 
saturation at the third codon position (cdp) for all PTLVs similar to 
that observed by others. "'"^^ However, substitution saturation was not 
observed in the first and second cdps and these sites were thus suitable 
for estimating posterior evolutionary rates and divergence dates of 
PTLV by using Bayesian analysis with an Markov chain Monte 
Carlo algorithm. Phylogenetic analysis of concatenated gag-pol-env 
first and second cdp sequences from PTLV-4 and prototypical PTLVs 
showed strong clustering of STLV-4 with HTLV-4 ( 1863LE) and sepa- 
rate lineages for PTLV-1, PTLV-2 and PTLV-3 (Figure 4C) consistent 



with analysis of the shorter gene regions. To estimate the possible date of 
introduction of STLV-4 into humans, we inferred divergence dates 
using a Bayesian framework and four combinations of clock models 
and tree priors. Bayes factor analysis showed that the relaxed clock 
model with a Yule process speciation model fit the data better for the 
TMRCA estimates than the relaxed clock model with a constant tree 
prior and both molecular clock models (21nBF=97.23). Using these 
methods, we found that the PTLV posterior mean evolutionary rates 
ranged from 4.68X 10"^ to 13.13X 10"^ substitutions/site/year (Table 3), 
which are consistent with previous estimates.*''^''^ Divergence of the 
progenitor of the PTLV-4 and PTLV-2 lineages using the relaxed clock 
model and a Yule tree prior occurred about 131 000 ya and were some- 
what higher using the three other models (Table 3), but were within the 
TMRCA date ranges previously determined by our group using analyses 
of separate gene regions.'^ In contrast, TMRCA estimates for the PTLV- 
4 ancestor are relatively more recent occurring about 3300-8400 ya, 
suggesting a relatively recent introduction of STLV-4 into humans 
(Table 3). 

STLV-1 and STLV-4 co-infection in one wild gorilla 

Sequence analysis of the partial env fragment (326-bp) obtained using 
PTLV generic PCR primers from goriUa 9584 showed that this animal 
was infected with both STLV-1 and STLV-4. The env sequences of this 
goriUa shared 98% nucleotide identity to STLV-1 from olive baboons 
{Papio anubis), which have a broad geographic distribution, including 
Cameroon. The 9584Ggo sequence also showed 99% nucleotide iden- 
tity to HTLV-1 subtype F env sequences from Gabon,^' which was 
confirmed by phylogenetic analysis (Figure 6). We consistently 
obtained STLV- 1 env sequences in repeat PCR reactions using dried 
blood spot DNA from this goriUa. However, examination of all 
sequences obtained from this animal showed no evidence of hetero- 
geneity as would be expected in a mixed infection, suggesting that all 
primers used for testing samples from this goriUa may have partial 
selective specificity to either STLV-1 or STLV-4, which may also 
explain why we were unable to amplify aU STLV-4 sequences from 
this animal using HTLV-4-specific primers. 

Sequence accession numbers 

AU new STLV-4 and G. gorilla mtDNA sequences generated in the 
current study have been deposited at GenBank with the accession 
numbers (KF029402-KF029433). 

DISCUSSION 

Prior to this report, HTLV-4 was the only known human retrovirus 
without a recognized animal reservoir. Previous efforts to identif)? 
HTLV-4-like viruses in a variety of African NHPs''"^' may have been 
unsuccessful because they relied on HTLV-1 -based serological assays 
and generic PTLV PCR assays with presumably limited sensitivities for 
STLV-4 detection. In addition, previous studies were limited in that 
they largely examined monkeys rather than gorillas or other apes for 
PTLV infection, since blood specimens are not as readily available 
from ape species.""^' The predominance of ape specimens from pre- 
vious simian retrovirus studies are non-invasively collected feces and 
urine samples^"'^'^^ and while SIV and SFV have been successfully 
detected in these body fluids, the detection of PTLV in these samples 
has not yet been reported.^"^^'^' 

A combination of screening data, geographical distribution, beha- 
vioral data and phylogenetic analyses provides strong evidence that 
gorillas are the natural host of STLV-4 and the reservoir for the known 
human case of HTLV-4. Through the largest NHP screening to date 
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Figure 4 Identification of STLV-4 in gorillas from Cameroon. Phylogenetic trees inferred using Bayesian analysis in the BEAST program and a relaxed molecular clock. 
(A) Topology inferred using a 225-bp alignment of tax sequences; (B) topology inferred using a 376-bp alignment of first and second coding position po/ sequences; (C) 
topology and divergence dates inferred using first and second codon positions of concatenated gag, po/ and eni/ sequences. Scale bars show units of time and for the 
concatenated sequences is years to the most recent common ancestors backwards in time from PTLV sequences. Posterior probabilities > 0.7 are provided at branch 
nodes. 
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Figure 5 Detection of nucleotide substitutions and APOBEC-induced hypermu- 
tation in alignments of complete and partial PTLV genomes. Analysis performed 
using tlie Hypermut and Highlighter programs at the HIV sequence database at 
the Los Alamos website (http://www.hiv.lanl.gov/). Sequence differences and 
evidence of hypermutation were determined by comparison with the complete 
genome of (A) HTLV-4{Caml863LE) or (B) STLV-3(CT0604) and are highlighted 
as follows: A, Green; T, red; G, orange; C, light blue; lUPAC, dark blue; gaps, gray; 
circle, APOBEC; diamond G>A. Base number is shown on x-axis, y-axis has no 
units and is used for only locating reference and query sequences. 



for evidence of PTLV-4, we detected STLV-4 in three wild and three 
captive gorillas, and found no evidence of PTLV-4 from a large series 
of other samples from a wide range of NHP species. While our know- 
ledge of the exact geographical origins of the three captive infected 
gorillas remains imperfect, when combined with the distribution of 
the infected wild gorillas, the data demonstrate that the individual 
infected with HTLV-4 lived in an endemic region of gorilla STLV-4. 
Furthermore, the HTLV-4-infected individual reported contact with 
gorillas through hunting, a behavior reported by only 11% of male 
respondents.^"* Given his self-reported behavior and the detection of 
PTLV-4 in gorillas from the area where the HTLV-4-infected indi- 
vidual lived and hunted, it is likely that his infection was the result of 
direct contact with gorillas through hunting. Phylogenetic analyses, 
including results from complete STLV-4 genomes from three gorillas, 
show that the gorilla STLV-4 sequences are nearly identical to that of 
the HTLV-4 case, lending further support to the conclusions that 
STLV-4 is endemic in gorillas, and that the HTLV-4 infection repre- 
sents the cross-species transmission of STLV-4 from gorilla to human. 
While the possibility that STLV-4 in gorillas resulted from the cross- 
species transmission from a third species cannot be ruled out, it is 
unlikely given the range of NHP species screened in Cameroon. 

Our results put to rest the debate over ancient origins of HTLV-4 
in humans. While the PTLV-4 lineage is older than the other HTLV 
lineages, the data indicate that human infection resulted from a 
relatively recent cross-species transmission of a gorilla virus, rather 
than an ancient infection in human ancestors. Our findings support 
the importance of gorillas in the emergence and evolution of retro- 
viruses and other human infections. In addition to HTLV-4, con- 
tact with gorillas is now thought to be associated with the origin of 
human infections with HIV-1 group P^'^^ and with sporadic human 
infections with STLV-lggo'^ and SFVgor.^'^'^''''''^ Interestingly, all 
variants of these novel retroviruses emerged in Cameroon where a 
long history of contact with NHPs and hunting of gorillas exists. 
Furthermore, recent work on the origins of malignant malaria indi- 
cated that gorilla Plasmodium was the source of human Plasmodium 
falciparum.^^ It now appears evident that gorillas should be con- 
sidered keystone reservoirs for human disease emergence along with 
other taxa such as chimpanzees and some groups of bats. Whether 
this is because of their phylogenetic proximity to humans or other 
elements of their population genetics, ecology or behavior remains 
to be determined. 

Our results also have conservation and clinical implications for 
gorillas, an increasingly rare close relative of humans. Gorillas are 
endangered and protected by law throughout their ranges, and the 



Table 3 PTLV evolutionary time-scale inferred using Bayesian analysis and first and second codon positions of concatenated gag, po/and env 
genes 



Clade/parameter 


Molecular clock (Yule speciation) 


Molecular clock (constant size) 


Relaxed clock (Yule speciation) 


Relaxed clock (constant size) 


PTLV root 


445156(352269-570416) 


466008 (367 842-591 500) 


233 060 (109 050-388 930) 


488 490 (234 220-853 200) 


STLV-5(MarB43VPTLV-l 


127368(100504-162 740) 


131065 (101099-165651) 


97953 (51554-156470) 


127921 (51561-233570) 


PTLV-1 


46449 (40000-57408) 


47 503 (40008-57871) 


48211 (40004-58209) 


49300 (40001-58525) 


HTLV-l(Mel)/PTLVla,b 


37886(30103-48314) 


38571 (30435^8776) 


35936(19533-51100) 


35754 (21608-51840) 


PTLV-4/PTLV2 


236731 (185905-303 707) 


244886(191135-311915) 


131072 (53805-233 563) 


215131 (94906-388425) 


PTLV-4 


3363 (1642-5604) 


3381 (1668-5712) 


3859 (749-8434) 


3707 (960-7495) 


PTLV-2 


145316(114368-187 615) 


150804(115280-190392) 


85261 (33005-147 580) 


122100 (47 534-217050) 


STLV-2 


46853 (34726-61419) 


47945 (35762-63054) 


27641 (5 592-53 657) 


33877(10222-53909) 


HTLV-2 


22 731 (17279-29336) 


23022 (17451-29 565) 


25775(15005^2 205) 


27 195(15050-42007) 


PTLV-3 


124744(97491-162517) 


128519 (98309-164636) 


102 100(43558-173160) 


138500(53454-247130) 


Mean substitution rate 


5.14X10"^ (3.94-6.3 XlO"^) 


5.02X10"^ (3.85-6.19 XlO"') 


8.71 X 10"' (4.68-13. 13 X 10"') 


5.314X10"' (2.85-7.99 XlO"') 


a-parameter (r-distribution) 


0.671 (0.358-1.014) 


0.669 (0.345-1.03) 


0.668(0.351-1.013) 


0.668 (0.352-0.984) 
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Figure 6 STLV-1 co-infectionof a wild gorilla from Cameroon. Identification of STLV-1 genotype in a wild gorilla using maximum likelihood phylogenetic analysis of 325- 
bp envelope sequences. Reference sequences were available at GenBank and accession numbers are provided in sequence names. One thousand non-parametric 
bootstrap replicates were used to assess the strength of the branching order. Bootstrap values ^70 are provided at branch nodes. Scale bar is substitutions per site. 
PTLV-1 subtypes are provided. Country codes are: Gab, Gabon; Jap, Japan; FrGui, French Guinea; FR, France; Can, Canada; Germ, Germany; Mel, Melanesia. Simian 
codes are: sm, sooty mangabey; bab, baboon; Ggo, gorilla; Tan, tantalus monkey. 



opportunities for synergetic messages from health, conservation and 
law enforcement authorities should be encouraged to ensure the max- 
imum benefit from any communication is obtained and that messages 
are harmonious. Considering that people continue to hunt apes in 
central Africa, greater efforts are needed to educate hunters on the 
risks of zoonotic infections to prevent the emergence of new microbes. 
The role that gorillas play in the emergence of human infectious 
disease provides an additional justification for prioritizing efforts to 
minimize hunting through strong enforcement of wildlife hunting 
regulations and work to develop appropriate and accessible animal 
proteins in regions where gorillas are endemic. 



While the clinical outcomes of STLV-4 infection in gorillas remain 
preliminary, the findings suggest that STLV-4 may be transmitted 
between gorillas in captivity and show that STLV-4 may lead to 
chronic infections in gorillas. While the significance remains unclear 
for now, the evidence of paralysis and neurological symptoms in 
Nyum are intriguing because of the known potential for similar symp- 
toms associated with HTLV- 1 and HTLV-2 in humans. Nonetheless, it 
is now clear that SIV can lead to serious disease in chimpanzees fol- 
lowing chronic infection,^ and leukemia and lymphoma have been 
reported in baboons infected with STLV-1, so the potential for 
STLV-4 to also cause disease cannot be excluded. The high pVLs 



Emerging Microbes and Infections 



A gorilla reservoir for HTLV-4 

M LeBreton etal 



11 



detected in some gorillas, including Nyum, may help explain the ease 
of transmissibUity of STLV-4 and the disease potential in these ani- 
mals. However, the two captive gorillas with high pVLs do not appear 
to have any clinical evidence of leukemia or lymphoma or other 
HTLV-l-like disease symptoms. It is unclear if the high pVLs in gorilla 
Nyum are associated with the possible neurological symptoms reported 
in this animal and wiU require testing of additional animals. In contrast, 
pVLs in the HTLV-4-infected hunter were comparatively much lower 
than those in the gorillas. It will be important to test close contacts of 
this hunter to determine if these low pVLs are associated with person- 
to-person transmission of HTLV-4. 

Our finding of a broad range of Western blot profiles in the STLV-4- 
infected gorillas using HTLV-1/2 antigens suggests that these current 
serological assays are likely missing or misclassifying PTLV-4 infec- 
tions in people and NHPs. For example, the HTLV-4 Western blot 
profile and those for the majority of the STLV-4-infected gorilla speci- 
mens from Shai were HTLV-2-like, suggesting that people with this 
Western blot profile, without follow-up PGR and sequence analysis, 
could be infected with PTLV-4 and not HTLV-2. While adequate 
specimen volumes were not available to evaluate detection of STLV- 
4 using available HTLV-1/2 EIAs, these Western blot results indicate 
that better serological assays are clearly needed to address this defi- 
ciency. However, the sensitivity and specificity of the HTLV-1/2 
Western blot assay for detection of PTLV-4 infection is currently 
unknown.'" The use of HTLV-4-specific screening PGR assays may 
have Umited the detection of PTLV-4-related infection in some ani- 
mals whose STLV-4 sequence is divergent in the qPGR primer regions. 
The development and use of the highly degenerate PTLV qPGR assays 
combined with the development of PTLV-4-specific, or improved 
serologic assays for broader PTLV detection, will all aid in better 
understanding the geographic distribution of PTLV-4 infection in 
humans and gorillas. Gombined with our finding of a wide distri- 
bution of STLV-4 in gorillas across the forests of Gameroon, the 
absence of HTLV screening in Gameroon and the use of non- 
HTLV-4 antigens for serologic detection, suggests a potentially 
broader human exposure to STLV-4 in Gameroon and that additional 
HTLV-4-infected humans may yet be identified. 

We also identified dual STLV-1 and STLV-4 infection in one wild 
gorilla (Ggo9584). Others have reported dual STLV-1 and STLV-3 
infected monkeys in Gameroon,''" suggesting that although rare, mul- 
tiple infections can occur in this region of high viral and primate 
diversity. The significance of these co-infections is currently unknown, 
but could lead to genetic recombination as occurs in other dual retro- 
viral infections. The STLV-1 in this gorilla clustered with HTLV-1 
subtype F from persons Uving in Gabon''^ and thus, could be the pri- 
mate origin of these human infections since the geographic range of G. 
gorillu extends west from the Gongo River into Gameroon, Gabon, 
Gentral African RepubUc and mainland Equatorial Guinea. Gorilla 
9584 was from near the southern border of Gameroon, a border which 
Gameroon also shares with Gabon and primate bushmeat hunting and 
consumption are known to occur in Gabon, further supporting this 
possibility. Thus, expanded testing for STLV-4 and HTLV-4 in Gabon 
and other African countries where gorillas exist may be warranted. 

The recent discovery of two new retrovirus species in humans 
(HTLV-3 and -4) and the ongoing detection of new cases of cross- 
species transmission of retroviruses from various animal hosts,'"^"'""'^ 
indicates the regularity with which humans are acquiring novel retro- 
virus diversity. Given the devastating nature of the global HIV pan- 
demic and the prior emergence of two species of HTLV that have 
spread globally (HTLV-1 and -2), the importance of monitoring for 



novel retrovirus outbreaks in human populations is paramount. 
Determining the animal reservoirs of zoonotic retroviruses, monitor- 
ing cross-species transmission of these viruses and others, and follow- 
ing individuals with primary infections all remain primary tools for 
the identification of new viruses with pandemic potential. Work to 
monitor individuals exposed to gorillas through hunting and other 
activities, appropriate efforts to decrease such activities and ongoing 
laboratory studies to further characterize STLV-4 provide a clear path 
towards preventing HTLV-4 emergence as a human disease. 
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